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ABSTRACT 

The effect of atmospheric turbulence on the power fluctuations of large wind turbines 
is studied. The significance of spatial non-uniformities of the wind is emphasized. The 
turbulent wind with correlation in time and space is simulated on the computer by 
Shinozukas method. The wind turbulence is modelled according to the Davenport spectrum 
with an exponential spatial correlation function. The rotor aerodynamics is modelled 
by simple blade element theory. Comparison of the spectrum of power output signal 
between 1-D and 3-D turbulence, shows the significant power fluctuations centered 
around the blade passage frequency. 


INTRODUCTION 


One of the con s i derations associated with the 
integration of large wind turbines into an electri- 
cal network is the fluctuation of output power/ 
torque due to turbulence in the atmospheric wind. 
Information on the statistics of shaft power/torque 
fluctuations Is essential for the design of control 
systems for large wind turbines. 

The effect of one dimensional turbulence (i.e. 
only temporal variations about a uniform mean wind 
over the entire rotor plane) has been studied by 
Frost [Ref. 1] and der Kinderen et al [Ref. 2]. 

The one dimensional analysis shows how the high 
frequency oscillations 1n power due to the wind 
fluctuations can be attenuated by large rotor in- 
ertia. _ : 

In this paper we emphasize the significance 
of spatial turbulence on large diameter rotors. 

[See Fig. L] Observing the wind energy spectrum 
of turbulence in the atmosphere [Ref. 26] we see 
that the predominant energy is at the low frequency 
end. This low frequency is associated with large 
size eddies on the order of hundreds of feet. When 
the rotor diameter approaches the size of these 
eddies, then the effect of the blade chopping 
through this spatial turbulence is to produce fluc- 
tuations in the power output at the rotor frequency 
and its multiples. The study is subdivided into 
two distinct tasks. 

(i) a suitable model to describe the turbu- 
lence in the atmosphere and the simu- 
lation of wind velocities conforming 
to this model . 

(1i) a suitable rotor model to produce the 
fluctuating power output signal with 
the above wind input. 

WIND MODEL 


The objective of the wind model is to simulate 
the turbulent wind components. The model takes 
into account the structure of the atmospheric 
turbulence, and generates a turbulent wind signal 
appropriate to the terrain, the scale of turbu- 
lence' and the mean wind velocity. 


A comprehensive survey of the subject of atmos- 
pheric turbulence [Ref. 16,17,18,19,24,25,26,27] 
gives a good understanding of the mechanics of 
turbulence in the atmosphere and also empirical 
equations to describe the structure of turbulence 
close to the ground [13,14]. 

Later we will show that the horizontal compo- 
nent of wind turbulence makes the largest contri- 
bution to the output power fluctuations. Hence 
our attention is focused on the modelling of the 
horizontal component alone. 


The Davenport Model 


According to Davenport [14] the spectral den- 
sity function for horizontal component is given 


by: 
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where: x = 4000. n 

n = frequency (Hz) 2 2 
S(n) = Spectral density (m /s per Hz) 

= Mean Wind velocity at reference height 
of 10m(33‘) above ground (m/s) 

Cj = Ground roughness factor (akin to rough- 
ness coeff. for pipe flow). Suggested 
values: 0.0005 for open sea, 

0.050 for Urban area 


The shape of this function [Fig. 2] conforms 
to observed general results of experimental measure- 
ments [Ref. 14]. The value of the constant Cj 
can be adjusted according to the terrain to give 
a model of the atmospheric turbulence at a parti- 
cular site. 

Davenports analysis of recorded data also 
shows that an exponential correlation between 
points in space agrees well with experimental 
results. The Spatial Correlation function is 
given by: 
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Cor. = exp(-^^-^) 

where: V = Mean Wind velocity 

C = correlation constant 


( 2 ) 


Shinozuka has shown that with as low as 50 fre- 
quencies one can obtain a fairly good represen- 
tation of the spectrum. [Ref. 22 J This proce- 
dure can be expanded to 3 dimensions as follows: 
[Ref. 21,22,23] 


This function is plotted in Figure 3 for C= 6 . 
Taking the Fourier transform of this function 
with respect to we get a spectral description 
in terms of the wave number « 

k = ^ (^ads/m ) 

we get therefore, for spatial correlation: 


S{k) = gjui/ (3) 

Tr(a^a)^+k^) 

where: a = C (per m ) 

2tt V 


For a complete 3-dimensional simulation of 
the horizontal wind turbulence as a function of 
time and two spatial coordinates Y and Z in the 
plane of the rotor we use the composite spectral 
energy density function non-dimensional i zed by 
mean wind velocity 


F(t,y,z) 
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To obtain a good representation of the spectrum 
in the turbulent signal, the values of N, 
should be large. Direct evaluation of tne above 
function turns out to be extremely time consum- 
ing computation. 

By suitably rewriting the above equation the 
Fast Fourier Transform technique can be applied 
in place of the above triple summation. [Ref, 

23] 

For the present simulation the constants in 
the spectral density function are: 


Mean Wind; Vj^ ^ 3 ^ 32kmph(20 mph) 

Terrain Roughness Factor = = 0.005 


S(n,k^,k^) 

’^M,33 


4-C.p. 
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Where are the wave nunt>ers corresponding 

to the s,eparatie?i lengths in the y and z directions 
respecti vely . a and 3 are related to the correla- 
tion constant C. (Note: C could be different 
for correlation in y and z directions). 


Simulation of Wind Velocities 


A very comprehensive method of simulating 
a multi-dimensional multivariate process is des- 
cribed by Shinozuka [21]. A random process is 
simulated by a series of cosine waves at almost 
evenly spaced frequencies and with amplitudes 
weighted according to the spectral energy at the 
wave frequency. 

For the simple case of a random function of 
time with spectral energy function S{u^) the pro- 
cess is simulated as: 


F(t) = [2,S(w .)Au)]^cos(w.t+e ,) (5) 

L 3 J J 

j=l 

where: Aw = ^u 

rr 

w = Upper cutoff frequency, above which 
S(w) is practically zero, 
w. = jAw + 6w,6w is randomly small quantity 
<Aw/ 20 is introduced to avoid a peri- 
odic repetition of the frequency w. 

0 = Random phase angle 0<e<2Tr ^ 


Spatial Correlation 

Constants a=3= 0.098 per m(,03 per ft) 

Cutoff Frequency n^= 0.5 Hz 

Cutoff Wave Number = k = 2 rads/m(0,628 rads/ 
ft) (chosen to give a 
grid spacing of 1.5m(5ft) 

After the simulation the turbulent wind data 
is stored as tabulated information. The table 
consists of 512 time planes at 1 sec intervals. 

In each time plane the wind velocities are tabu- 
lated at grid points 1.5m(5ft) apart in both y 
and z directions and covering a total area of 
98x98m(320 x320 ft) in the y-z plane. (64 points 
in each direction) 

Instead of generating the turbulent signal 
as a function of t,y,z by evaluating the triple 
series a table look up is done in the tabulated 
wind data with simple linear interpolation. 

In the present work we evaluate the turbulent 
wind component at different radial stations on 
the wind turbine blade. Interpolation has been 
restricted to between time planes only. For tur- 
bulence at a particular spatial position we choose 
the value at the grid point closest to this posi- 
tion. 

Details of the wind model and computational 
techniques are given in Ref. 30. 

An example of the accuracy of the simulation 
method is shown in Fig. 4 where a comparison be- 
tween the Davenport Spectrum and a spectrum ob- 
tained by taking the Fourier Transform of a one- 
dimensional turbulent signal generated using Eq. 

5, shows excellent agreement. 

AERODYNAMIC MODEL OF ROTOR 


All of the theories available for analyzing 
propellers are applicable for analysis of windmill 
rotors [3,4,12,13,14]. Based on the methods 
developed for propellers, the performance of wind- 
mills can be analysed [5,6,8,9,10,11]. Specific 
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n^thod for analyzing the unsteady aerodynamics 
has also been developed [8,9], The study by Barlow 
[3,4] presents a fairly comprenensive understand- 
1ng of the forces on a propeller operating in 
a turbulent atmosphere. As pointed out earlier, 
this analysis can be equally adapted to the wind- 
mill rotor. 


Simple Blade Element Theory 

For the purpose of the present analysis a 
relatively simple model of the rotor is adequate. 
The Blade -element theory is adopted because of 
its simplicity and ease of computation. Instead 
of computing tip loss corrections, and induced 
velocities, we simply introduced an overall in- 
duced velocity factor. 

The resulting expression for the power de- 
veloped is very simple. Comparison of the results 
from this simple model with that of the Wilson- 
Llssaman program shows that the trend of the 
results is the same. By a suitable choice of 
the induced velocity factor one can attain a speci- 
fic operating point defined by the tip speed ratio 
A and the power coefficient C , for a particular 
pitch setting angle e . ^ 

Although the mode? is very simple, it still 
retains all the nonlinear features of the wind 
input versus the power output. 

In analyzing the effect of wind fluctuation 
on output power we are interested more in the 
fluctuating component of power, than in the exact 
reproduction of the actual power. If one is inter- 
ested in the exact power at different operating 
points the more exact methods detailed above can 
be used. 


Integration of the Wind Model & Rotor Model 
The aerod ynamic model of the rotor consists 


of calculating the forces and moments on each 
element of the rotor from the resultant velocity 
at that station. The free stream wind is modeled 
as the sum of the mean wind velocity including 
the effects of wind shear and the turbulence com- 
ponent as generated by the wind mode!. That is 


V = V (z) + v'(t,y,z) (7) 

In terms of the freestream turbulent wind com- 
ponents u',v' and w' the fluctuating power 
coefficient is shown to be (Ref. 30) 

0 ^ (i+X^x^) 


{^'sine - w’coso + uj_ Ax]xdx (S) 

U U U 

0 = Qt = angular position of rotor 

A = = tip speed ratio 

U 

From this equation we see the dominance of the 
uj term which is multiplied by A the tip speed 

ratio. 


In the present model we only use the horizontal _ 
component of wind turbulence and treating the ^ 
free stream wind velocity to be only in the axial 
direction we can easily find the resultant fluctu- 
ating wind input at each point along the rotor. 

As the rotor goes around in space each point on 
the rotor is at different spatial points at 
different times. Knowing the spatial point coor- 
dinates and time, the wind model can generate 
the turbulent velocity which then is added to 
the mean wind velocity to generate the result- 
ing local wind velocity at that rotor station. 

It can be observed, that the output power 
would be fluctuating even if the input wind has 
a frozen turbulence. That is, there is correla- 
tion with respect to space only. An analysis with 
frozen turbulence for a simple element of the 
rotor is presented in [Ref. 30]. 

RESULTS AND DISCUSSIONS 


The aim of the present study is to show the 
important contribution of the spatial non uni- 
formities of atmospheric turbulence on the output 
power fluctuations of a large wind turbine. 

Using the Davenport spectrum and the exponen- 
tial correlation function, with Shinozukas method 
of generating a random signal a complete table 
of turbulent wind data haA been generated. 

The next step is to input this wind model 
to the rotor model and analyze the resulting power 
output. In order to emphasize our point about 
the contribution from the spatial turbulence 
effects, we first study the simple one dimensional 
wind turbulence Input. 

In Fig. 5 is shown the unsteady power output 
with a one dimensional wind input. It must be 
pointed out here that the power output obtained 
is an ideal one. This is the direct power 
obtained from the aerodynamic forces actingon 
the blade. The inert la of the rotor and the 
flexibility of the shaft have not been considered. 
For such an ideal case obserye the exact repro- 
duction of the wind fluctuations In the power 
fluctuation. 

Spectral analysis of the input one dimensional 
wind turbulence and' the output fluctuating power 
is shown in Fig. 6. As expected from our obser-^ 
vat ion of the time signals the spectra are identi- 
cal in shape. Comparison of the Davenport spectrum 
of Fig. 3 with the spectrum of the wind turbulence 
Fig. 6 shows that the envelope of the peaks follows 
the Davenport spectrum. Although the simulated 
spectrum shows the correct shape, due to the finite 
length of the generated time signal the spectrum 
is not smooth at the low frequencies. This could 
be remedied by generating a longer time signal. 

To show the effect of spatial turbulence we 
considered three different diameter rotors 30m, 

60m, and 75m (100 ft, 200 ft, and 300 ft). 

The results for the three different sizes 
of rotors are shown in Figures 7 to 9.* For com- 
parison in each case is shown the power output 
with a one dimensional wind input and above it 
the power signal with a three dimensional wind 
input. (By three dimensional we mean the turbu- 
lent wind component as a function of time t, and 
spatial coordinates y and z in the planes of the 
rotor). 


* All cases are 
same tip spee 


for 2-bladed wind machines at the 
d ratio 

U 
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Observe the high frequency fluctuations caused 
by the rotor chopping through the spatially non- 
uniform wind. The power spectrum of the time 
signals adequately show the important contributions 
from th^ spatial turbulence effects. Comparing 
the spectrum of the three different size rotors, 
observe how with increasing rotor diameter, the 
fluctuations caused by spatial turbulence dominate 
over the original Davenport spectrum due to the 
temporal fluctuations of the wind. 

The blade passage is also indicated in the 
figure to show that the spatial turbulence effects 
are concentrated around the blade passage fre- 
quency. 

Conclusion 

The study conducted in the present simulation 
has shown the importance of spatial turbulence 
on performance of large wind turbines. This aspect 
of the fluctuations is an important consideration 
in the design of pitch control systems. 

The above study can^ be easily extended to 
incorporate a more sophisticated rotor model, 
like the Wil son-Lissaman program or the vortex 
lattice method. In calculating the power output 
at the shaft the flexibility of the shaft and 
the rotor inertia can be introduced by writing 
the dynamic equilibrium equation for the entire 
rotor system In the presence of the applied aero- 
dynamic torque and the torque load on the shaft. 
Solution of this dynamic equation will then pro- 
vide the necessary shaft torque or power as a 
function of time. This extended model can be 
used to study the influence of rotor inertia and 
the flexibility of the shaft. 

This complete dynamic model can then be merged 
into an overall simulation of a wind turbine power 
generation system as is done in Ref. 23 and 29. 
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FIG. 1 THE PHYSICAL PROBLEM 



LOC( FREQ<H*) > 


FIG. a DAVENPORT MODEL FOR SPECTRUM OF 
HORIZONTAL COMPONENT OF WIND TURBULENCE 



REDUCED FREQ <n*l/V) 

FIG, 5 CORREUATION FUNCTION FOR SPATIAL 
CORRELATION OF WIND TURBULENCE 
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QUESTIONS AND ANSWERS 


J.P. Sullivan 


From: W*N. Sullivan 

Q: Please comment on the effect of very soft drive trains on the output power spectra. 

A: A simulation of this problem is planned for the near future at Purdue, 

From; Walter Frost 

Q: How is the coherence between frequencies of the wind simulated? 

A: With the exponential correlation function exp [- ] a - 6, 

From: F,W, Perkins 

Q: Why does the turbulence seem to drift in a preferred direction, low right to high 

left? 

A: I don*t know^ hut it may he just a figment of color image processing. 
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